隐式Newmark法分析波动问题精度的探讨

本文研究了用平均加速度的Newmark积分格式求解一维平面弹性波动问题的精度。由于平均加速度的Newmark积分格式没有振幅衰减、不产生数值阻尼,因此采用此法分析波动问题会产生较强的高频振荡,尤其是加速度的计算值误差甚大。本文对防护工程中常见的平台型荷载和升压三角形荷载作用下的波动问题,应用有限元分析的实例表明:采用有数值阻尼的Newmark积分格式是改进有限元求解波动问题的一个有效方法,可以大大地提高应力和加速度求解精度。文中推荐的有数值阻尼Newmark积分格式的积分参数(=0.9,=0.49)可供工程中土与结构动力分析时参考使用。     

推进－迭代算法计算效率的研究

本文对求解三维定常超音速动性流场的一次空间推进，在每一个推进站沿伪时间层局部迭代的推进－迭代算法作了进一步的研究．在每一推进站（侧向平面）沿伪时间层局部迭代时，给出了四种不同的隐式迭代方法，即沿侧面两个方向（法向和周向）全用隐式；法向隐式而周向采用Ｇａｕｓｓ－Ｓｉｌｄｌｅ来回扫描迭代；法向隐式而周向显式及以系数矩阵谱半径代替系数矩阵的简化标量隐式算法．用这四种算法模拟了三维球锥黏性绕流，给出了四种不同算法的计算效率和收敛特性比较．     

无条件稳定数值积分方法在拟动力实验中的应用研究

用无条件稳定的隐式α－方法代替目前的条件稳定的显式方法，通过数值迭代控制方法实现了拟动力实验．在此基础上进一步应用子结构技术，解决了刚度大、自由度较多结构采用传统的条件稳定数值积分方法进行拟动力实验的问题．应用不同数值积分方法进行的三层模型结构的实验结果表明，无条件稳定的数值积分方法的实验精度优于条件稳定的数值积分方法的精度．文中还就实验误差的控制补偿问题进行了讨论，提出了采用瞬时割线刚度进行补偿的方法     

On numerical solution of dynamic problems in the theory of reinforced shells

Dynamic problems for cylindrical shells reinforced with discrete ribs are examined. A numerical algorithm based on Richardson extrapolation is developed. Specific problems are solved, and the results are analyzed __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 5, pp. 50–56, May 2006.     

一种用于分散化滤波的融合算法及其在组合导航系统中的应用

提出了一种局部估计相关时的最优融合算法，其特点是加权矩阵是对角阵，即全局估计中的各分量分别是局部估计中相应分量的线性组合，在此融合算法中局部估计是最优的，并且不存在由主滤波器向子滤波器的信息反馈，因此具有较好的客措性。给出了最优加权矩阵的存在性和唯一性的证明，并通过与现有类似融合算法的对比，表明该算法具有计算量小的优点。最后，以组合导航系统为例作了仿真计算，结果表明该融合算法是可行的。     

基于非结构化同位网格的SIMPLE算法

通过基于非结构化网格的有限体积法对二维稳态Navier—Stokes方程进行了数值求解。其中对流项采用延迟修正的二阶格式进行离散；扩散项的离散采用二阶中心差分格式；对于压力-速度耦合利用SIMPLE算法进行处理；计算节点的布置采用同位网格技术，界面流速通过动量插值确定。本文对方腔驱动流、倾斜腔驱动流和圆柱外部绕流问题进行了计算，讨论了非结构化同位网格有限体积法在实现SIMPLE算法时，迭代次数与欠松弛系数的关系、不同网格情况的收敛性、同结构化网格的对比以及流场尾迹结构。通过和以往结果比较可知，本文的方法是准确和可信的。     

结构拓扑优化设计的三角网格进化法

针对进化或拓扑优化方法的不足，提出了一种基于遗传算法的新型进化式拓扑优化方法－－三角网格进化法，该方法不仅能够同时进行拓扑，形状与截面变量优化设计，而且在优化过程中实现了退化和进化的统一，提高了优化效率。另外本文还首次对结构类型变量进行了优化计算，取得了有益的结果。最后几个数值算例证明了本方法的可行性和有效性。     

一种改进的分块隐式数值方法

提出了一种改进的分块隐式数值方法，在贴体坐标和交错网格下以逆变速度分量和压力Ｕ，Ｖ，Ｗ，ｐ为基本求解变量，由此克服了原分块隐式数值方法求解复杂边界流动时的困难．９０°弯管流动数值计算初步表明，本文提出的方法合理、可行     

On the solution of plane boundary-value problems of elasticity in a rectangle

An algorithm for solving plane boundary-value problems of elasticity for a rectangular domain is expounded. The algorithm is based on a complex-valued representation of the general solution to the differential equations of the plane problem and on the use of Lagrange polynomials to satisfy the boundary conditions. The algorithm can quite easily be implemented in a computer program. This is probably the simplest way of solving boundary-value problems of this class __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 1, pp. 97–102, January 2006.     

部分输入未知条件下结构动力复合反演的分解算法

在充分利用部分输入已确知而部分输入未知的激励特性的基础上，提出了结构动力复合反演的分解算法，该算法从源头上消除了迭代过程中参数识别与荷载反演的相互影响，降低了问题的计算规模。对于线性参数系统，该算法不经过任何迭代计算即可一次性完成结构参数识别及荷载反演。将其与松弛法结合，可解决非线性参数系统的识别问题，与文献[4]的方法比较，其收敛速度有显著提高。